Dynamic Compressibility and aging in Wigner crystals and quantum glasses 
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We study the non-equihbrium linear response of quantum elastic systems pinned by quenched 
disorder with Schwinger-Keldysh real-time techniques complemented by a mean-field variational 
approach. We find (i) a quasi-equilibrium regime in which the analytic continuation from the 
imaginary-time replica results holds provided the marginality condition is enforced; (ii) an aging 
regime. The conductivity and compressibility are computed. The latter is found to cross over from 
its dynamic to static value on a scale set by the waiting time after a quench, an effect which can be 
probed in experiments in e.g. Wigner glasses. 
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The dynamics and transport properties of glasses are 
the object of current theoretical and experimental inter- 
est III y, y{ . The slow approach to the static limit and 
the accompanying aging phenomena 1] observed exper- 
imentally in classical glasses is captured by a variety of 
models j3| • A natural explanation for the slow dynamics 
is the existence of a special organization of an exponen- 
tially large number of metastable states. 

Less is known when glasses evolve at very low tem- 
peratures and quantum fluctuations become irnportant 
as occurs in Wigner 0, 0| and Coulomb 0, S B' 
glasses, as well as in spin ^3 ^^'^ other systems 
at low temperature. Aging in the transport properties of 
Coulomb glasses was reported recently p| . Theoretically, 
electronic glasses have been studied with numerical simu- 
lations in a classical limit or with the imaginary-time 
methods combined with the replica trick applied 
to the isolated model. Strictly, the latter allows one to 
describe the statics and, via an analytical continuation 
to real time, the equilibrium dynamics with infinitesi- 
mal dissipative coupling to a bath. However, it is of 
a somehow broader use since it has also given access to 
dynamical quantities such as the finite frequency conduc- 
tivity |l2L Il3l Il4| , and instanton calculations allowed to 
relate, in the presence of a bath, the imaginary time so- 
lution to the ultra-slow dynamics (quantum creep) 
Still, a rigorous treatment of a full non equilibrium relax- 
ation requires special techniques devised to deal directly 
with the real-time dynamics of dissipative quantum sys- 
tems Even in the infinitesimal coupling limit 
these are needed to ascertain the validity of analytic con- 
tinuations to real-time, especially for glassy systems. 

The compressibility of electronic glasses has been mea- 
sured recently 18, Ujl, |2(J. Being of thermodynamic na- 
ture, in equilibrium it should only depend on the statics 
of the problem. Conversely, its time dependence reflects 
an out of equilibrium relaxation. Within an imaginary- 
time (Matsubara) variational calculation [12 . Il.l . l2 ij the 
static, zero frequency, compressibility of a disordered 



Wigner crystal was found to be non-zero and identical 
to the one of the pure system. In contrast, the analytic 
continuation yielded a vanishing result even in the small 
real frequency limit. This hinted to the fact that aging 
effects could be present M, 121*] though a firm conclusion 
was clearly beyond the imaginary-time calculation that 
assumes equilibrium at the outset. 

In this Letter we study the out of equilibrium relax- 
ation of a disordered Wigner crystal with the Schwinger- 
Keldysh (SK) technique [2II23. We find two two-time 
regimes; one in which the equilibrium result is recov- 
ered and another one in which the compressibility is re- 
duced and dominated by aging effects. Our calculation 
is performed within a mean-field-like variational approx- 
imation. We discuss its limits of validity as well as the 
relevance of our results for experimental systems. 

We model a disordered quantum crystal as 



H 



hU (x) c,„ , 



U{x)pqCos[Q{x - u{x))] , (1) 



where m is the mass of the particles, po is the aver- 
age density, Q = l-nja with a the inter-particle spac- 



ing and 



n and u are conjugate operators 



[u(a;), n(a;')] = ihS{x — x'). U{x) is a random poten- 
tial with Gaussian statistics, U{x) = and U{x)U{x') — 
A(a; — a;') with lS.{z) a function with finite range rf. This 
model describes a large class of systems including charge 
density waves (with a phase (j) — 27rM/a), Wigner 

crystals (upon generalization to a two cornponent vector 
u) 0, 0| , and Luttinger liquids in d = 1 • 

The compressibility k is defined as the change in den- 
sity in response to a change in chemical potential p. In 
linear response, it is given by the q ^ limit of the 
density-density correlator. For |^ the long wavelength 
part of the density is p(x) — pa ^ ~pa'S/u{x), and the 
equilibrium compressibility is k = lim^^o = 0) 



2 



with K{q,uj„ 



p^q Gc{q,oJn), the Matsubara fre- 



quencies, and Gc{q,ujn) 



where (...) and • • • denote thermal and disorder average, 
respectively. Within the replica variational approach |l2l | 



K.{q,LO„ 



Pmi^n + + Sl(l - 5n,o) + /(Wn) 



(2) 



where pm = rnp^^ Ei ~ Pm^'j, is a constant depending on 
disorder, ojp is the pinning frequency 0,113 and /(O) = 
0. K is independent of disorder and simply given by 



X = 1/T where T is the temperature (fcs = 1). How- 
ever, the response to a time dependent field is always 
zero, leading to x(a; ^ 0) = at variance with the static 
result. The difference is due to the existence of two de- 
generate ground states in the unperturbed system. In 
the static calculation one sums over both. The pertur- 
bation instead is unable to induce transitions leading to 
zero response. In a glass there may be no exact degen- 
eracy between the multitude of metastable states. Still, 
if they are close enough in energy and their coupling is 
sufhciently weak one expects a similar phenomenon. 



K = pI/c . 



(3) 



Alternatively, the real-time compressibility, i.e. the re- 
sponse to a time-dependent chemical potential, is given 
by the retarded linear response. Naively, this can be ob- 
tained from ^ by the standard analytical continuation 
iujn uj + iS that leads to [121 



K{q, oj) = plq^ [-Pm^'^ + cq^ + Si + /(t^)]" 



(4) 



where I{u) ^ 0) ^ 0. Note that even performing the 
limit LO ^ Q first while keeping q fixed one finds 



K = lim lim K{q,Lu) = , 



(5) 



in disagreement with the result in P)). 

The simplest example where such a difference arises is 
an isolated two level system {e.g. a spin). The static 
susceptibility in response to an external magnetic field is 



To analyze this issue, we study the real-time dynam- 
ics of model Q using the SK technique. Let us sketch 
the main steps of the calculation. The observables are 
computed using a path integral on two fields, u±, that 
live on the two sides of a closed time contour. The mea- 
sure is given by e^^'^ where Sk is the (dimensionless) 
SK action Sk = ji[S{u+) — S{u-)] and S{u) is the stan- 
dard action for a system described by To take into 
account non zero temperature and dissipation we couple 
the system to a thermal bath of independent harmonic 
oscillators 10 . By integrating them out we induce a cou- 
pling between the fields u±. The average over disorder 
is easily done without replicas and has a similar (cou- 
pling) effect. In terms of the fields u = ^{u+ + m_) and 
u = ji{u+ — U-) the resulting action, Sav = Sq + Sd, 
depends on R(z) = PgAg cos{Qz) with Aq the Fourier 
transform of the disorder correlation at wavevector Q, 
and reads 



^0 





-Pmuo'^ + cq^ — irjuj 



-Pm^ 

-77ftwcoth(|^ 



cq + ir\L>J 



^q^uj 



2h? 



ct=±1,ct'=±1 



(6) 



with /^^ = J J 1^. The terms proportional to i] 
arise from the coupling to the bath and represent Ohmic 
dissipation. For jj infinitesimal one recovers the intrin- 
sic dynamics of the quantum system. For the pure 
system (R = 0) this action yields the equilibrium re- 
sponse obtained by the analytic continuation of the Mat- 
subara representation. Sd can be rewritten as 5*^ — 

^^W^ Lw sin((3Ma;t/2) smiQhu^t' f^) cos[Q{ua:t-Uxt')]; 
it is clear that the quantum action crosses over to the dy- 
namic Martin-Siggia-Rose action when ft 0. 

The disorder-averaged correlation and linear response 
are defined as Cx-x'{t,t') ~ ^{'^tt'^x't' + ^xt^t't'' 
{uxtUx't') and Rx-x'{t,t') 



{uxtiux't'), with Rx-x'it,t') for t < t' and {uxtUx't') van- 
ishing because of causality. The brackets represent here 
an average with the weight e"'^"" . In the absence of dis- 
order, C and R are stationary, i.e. depend only on t — t', 
their Fourier expressions are 



{pmUj"^ — cq^Y + r\^iJ^ ' (7) 

2 I „„2 „■ ,1-1 



^q^'^) = {-Pra^ + Cq" - ITycj] 

and obey the fluctuation-dissipation theorem (FDT) 
2tanh(^) 



i?"(u;) 



ft (w — w' -f ie) 



(8) 
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In particular, Mmi?°(tj) = tanh(/3fiw/2)C°(w')- Since 
the long wavelength part of the density and the current 
are respectively p{x) — po = —Po'^u{x,t) and j{x) — 
PQdtu{x,t), in linear response the compressibility and 
conductivity are given by 

K(q,uj) = q^plRq{u}) , cr(w) = -iujpl Rq=Q{u) . (9) 

Inserting the free values in lO one recovers 10) for the 
static compressibility and <j{uj) = p^/iTj — iwpm) for the 
static conductivity which, in the limit 77 — > 0, reproduces 
the Drude form, Recr(tj) = ^S{uj). 

In presence of disorder, one derives self consistent 
mean-field equations for Cq{t, t') and Rq{t, t') where t — Q 
is the time when the system is quenched into the disor- 
dered state and set in contact with the bath. We show 
only the equation for the response (see i28| for details) 

{p„,d''t + rjdt + cq^)Rqit, t') = 5{t - t') 



dsi:{t,s)[Rq{t,t') - Rq{s,t')\ , (10) 



^lYaV'[B(t,s) + ihR(t,s)] 



with the self-energy E(i, s) = 

and V{z) = — ppAge^^'^^^. The tilde denotes an inte- 
gration over q [e.g. R{t,s) — J^Rq{t,s)] and Bq{t,t') = 

Cqit,t)+Cqit',t')~2Cq{t,t'). 

Analysis of this equation in the lon g tim e limit t, t' 
00 shows that, as in the classical case |29( and the quan- 
tum dissipative p-spin model [lil Il7| , the model exhibits 
two two-time regimes. First, for fixed t — t' the two- 
time functions are stationary, Rq{t,t') — > rq{t — t') and 
Bq{t,t') bq{t - t'), and the FDT © holds. For more 
separated times, there is an asymptotic aging solution, 
Rq{t,t') Rq{t,t'), as discussed below. The equation 
for rq{t~t') is obtained from 111) 1 1 bv a careful separation 
of time scales along the lines of 'SO^ . Its solution reads 

rq{uj) = [cq^ + M - irjuj - pmUJ^ - S(w)]"-^ , (11) 

with the self-energy I](w) = dT{e-"^^ ~ l)S(r) and 

E(t) J2 ^'^y'[HT)+ificr?{T)] . (12) 

The constant M = limt-,^ J* dsE(t, s) - J^'^ dT'T.{T') 
the so called anomaly, arises from the contribution of the 
aging time-scales to the FDT regime. Using FDT, 



m = / 2( 



l_e-'-t);i[l + 2/B(co)]Imf(^) , (13) 



where fsii^) — l/{e^^^'^ — 1). The last self-consistency 
condition follows from matching the FDT regime with 
the aging one. Taking t' t~ in (|10|l yields {cq^ + 
M)Rf{t,t~) = Y.^{t,t')rq{uj = 0), and the existence 







time 



FIG. 1: Sketch of the variation of the chemical potential; do 
(solid line) and ac (dashed line), for a system quenched at 
time t = Q. The density is measured at time t. 



of a non vanishing aging solution with R^{t,t') 
t' — > requires the 'marginality condition' 



1 = -4y"(5oo) / ( 

J a 



cq^ + Aiy 



as 



(14) 



where 600 — linit^ooK^)- O'^^ can explicitly check [23 
that the solution in the FDT regime coincides with the 
analytical continuation of the saddle-point solution of the 
replica variational approach to the Matsubara action per- 
formed in appendix D of More precisely, Gc{q,(^n) 
identifies (after analytic continuation) with rq{uj). Sim- 
ilarly, /(iw„ -> w + «(5) -> /(cj) = -E(w), El M, 
B — > 600 where B is defined in (31) in ,12|. Importantly, 
this correspondence holds only if the replica symmetry 
breaking scheme is the one using the marginality con- 
dition (|14|l . This choice, advocated in l^l for being the 
only one leading to a gapless conductivity, was often used 
since and is hereby fully justified within the SK formal- 
ism. Note that in d — 1 and for 77 = model has 
an additional stable 1 step replica symmetry breaking 
solution that is then discarded in the Matsubara treat- 
ment. The SK formalism also allows one to obtain the 
response for finite j]. The low frequency behavior of 
the conductivity changes from Re(T(a;) pquj'^ / (mujp) 
to Rea(w) ~ ^w3/2/u;3 ^ 

We now turn to the compressibility, i.e. the (linear) 
response of the density to a small change in chemical 
potential of amplitude dp applied between times tw and 
t (see Fig.^l. One finds 



K[q;t,tw) = 



6p 



2 2 



cq^ 



(15) 



As t — increases one distinguishes two regimes. For 
LUp^ t — tyj <^ the response is dominated by the 
'FDT regime' and it is stationary; moreover M{t,tw) = 
M and the compressibility vanishes [^(g; t, t^) when 
q — > 0] recovering If, instead, t — ^ tw, the re- 
sponse is dominated by the aging regime and M{t, tw) — 

M - 4 /|^(^ V"iz)X[z]dz, where X{z) is the FDT vi- 



olation ratio |lfil IHfij i as a function of z = B. For rj 
infinitesimal, X{z) coincides with X —^ u, z ^ B{u) as 
found in the replica solution with H14|l 12]. B has an ag- 
ing form and its detailed scaling depends on the model. 
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In d = 1,2 one finds M{t,t^) ~ MF{h{t)/h{t^)), where 
F and h are scaling functions. When times are very sepa- 
rated t — tyj iuij ^(i, tw) tends to zero and, from ((T3|l . 
the compressibihty is the constant ||2J). At intermediate 
time scales a mass is always present and the compressibil- 
ity depends on the wavevector q. For a finite size system 
of size L one can estimate that the compressibility crosses 
over from being essentially zero to the thermodynamic 
one when 



{LjLf> M{t,t^)/M 



(16) 



[Lc = yc/pm'^p is the Larkin pinning length] . The char- 
acteristic timescale to realize such evolution is the waiting 
time tyj. This potentially provides a direct experimental 
way to check for aging in these systems. Alternatively, 
one can apply an ac perturbation as typically done in 
experiments (see Fig. ^ and Q for a similar discussion 
on the ac magnetic susceptibility of classical spin-glasses) 
and obtain different results by tuning the period, r, of 
the perturbation. Short r's select the FDT contribution 
and thus a vanishing compressibility while very long r's 
do not erase the contribution from the aging regime yield- 
ing a constant compressibility as in Q . The crossover in 
r is typically of the order of . 

The above results are consistent with the picture that 
upon a change of chemical potential the system is first 
trapped in quasi-equilibrium in a metastable state. Since 
these states are pinned the charge cannot fluctuate and 
the compressibility is essentially zero. On the other 
hand if the change in chemical potential is maintained 
for a long time, the system explores other metastable 
states, the charge being allowed to change in the process, 
thereby leading to a finite response. The mean- field so- 
lution presumably overestimates the separation between 
metastable states. Transitions can occur through acti- 
vated processes due to quantum or thermal fluctuations 
(so called creep). Although treating such processes is 
difficult a possible modification of (|15|l is 

K{q; t, t^) = [qL{t^), h{t)/h{t^)] (17) 
c 

with F(0,?/) = (hF{x,oo) = F{oo,y) = 1. Classical 
creep arguments [S^ based on barriers growing as sug- 
gest L{t) = Lc + h(t) with h{t) = (Tln^^/e^ Extensions 
incorporating quantum effects, as in [l5|, are needed to 
complete this picture. 

In addition, the approach developed here allows one to 
generalize the mode-coupling theory to low-temperature 
glasses with quantum fluctuations and no quenched dis- 
order (see for a discussion of the classical Hmit). 
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